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THE DATA SINGULAR AND THE DATA ISOTROPIC LOCI FOR 

AFFINE CONES 

EMIL HOROBET 


Abstract. The generic number of critical points of the Euclidean distance 
function from a data point to a variety is called the Euclidean distance degree. 
The two special loci of the data points where the number of critical points is 
smaller then the ED degree are called the Euclidean Distance Data Singular 
Locus and the Euclidean Distance Data Isotropic Locus. In this article we 
present connections between these two special loci of an affine cone and its 
dual cone. 


1. Introduction 

Models in science are often expressed as real solution sets of systems of poly¬ 
nomial equations, namely real algebraic varieties. One of the most fundamental 
optimization problems that can be formulated on such sets is the following: given a 
real algebraic variety and given a general data point of the ambient space, minimize 
the Euclidean distance from the given data point to the variety. 

In order to solve this problem algebraically we examine the critical points of 
the squared Euclidean distance function. The number of such critical points ex¬ 
presses the algebraic degree of the complexity of writing the optimal solution to 
the distance minimization problem and it is called the Euclidean Distance Degree. 
This optimization problem arises in a wide range of applications, such as low rank 
approximations (Example]^, control theory (Example [T^, formation control (Ex- 
amplefl^, algebraic statistics (Example[T4|) and multiview geometry (Example [T^. 

Eor a general data point u the number of complex critical points is constant 
while, the number of real critical points is typically not constant for all general u. 
Eor example, if one of the critical points has a multiplicity, then the number of 
real critical points typically changes, this locus is called the ED-discriminant (or 
classically focal loci) and it was studied in P151151151IT^. 

In this article we want to discuss the locus (different from the ED discriminant) 
of exceptional data points u for which the number of complex critical points is 
smaller then the ED degree. There are three ways in which we can have different 
number of critical points than expected. The first reason is because a critical point 
may wander off into the singular locus of the variety. The study of this special locus 
was proposed by Bernd Sturmfels, first examples were developed by the authors of 
[3] and it was named ED data singular locus. In a similar fashion, the second 
case is when a critical point becomes isotropic with respect to the Euclidean inner 
product; this locus will be called ED data isotropic locus. In these two cases the 


2010 Mathematics Subject Classification. 14N10, 41A65, 55R80. 

Key words and phrases. Euclidean distance degree, Data Singular locus, Data Isotropic locus, 
ED Discriminant. 


1 



2 


EMIL HOROBET 


number of critical points is smaller then the ED degree. Finally, a data point can 
have infinitely many critical points, but this phenomenon is apparently recorded by 
the ED discriminant, so we do not deal with it in this article. A classical example 
would be that there are infinitely many critical rank 2 approximations of a matrix 
with two identical singular values. 

In this article we aim to describe the data singular and the data isotropic loci of 
affine cones. 

Acknowledgments. The author is grateful to Jan Draisma for his support and 
help and to Bernd Sturmfels for the theoretical insight. The author was supported 
by the NWO Free Competition grant Tensors of bounded rank. 

2. The special loci of data points 

In order to find the critical points algebraically, we consider A to be a variety in 
C" and we examine all complex critical points of the complexified distance, induced 
by the standard symmetric bilinear form, 

n 

du{x) = {u- x\u - x) = 

i^l 

with X G A™*, where A''®* denotes the locus of regular points of A, so we only 
allow those critical points that are non-singular. Let A C C" be an irreducible 
algebraic variety of codimension c with defining radical ideal I. If a: € A’’®® is a 
critical point of du, then the following holds: u — x T T^X. This later condition can 
be formulated as a: € A*'®® is a critical point of du if and only if all the (c-l-1) x (c-|-I) 
minors of the matrix 

/ u — x \ 

y dac^il) ) 

vanish, where JaCa;(/) is the Jacobian of / at the point x. 

We define the ED-correspondence to be the closure of the set of all pairs (u,x), 
such that x G A''®® is critical to du, and we denote it by fx C C” x C”. In other 
words Sx is the closure of the variety defined by 

rank(j“-,"^, )<c} 

We have two natural projection maps tti : Sx —> C" sending (u, x) to u and 
712 : fX —>■ C” sending {u,x) to x. Let SingA denote the singular locus of A, that 
is the set of all points of a; G A such that all the c x c minors of JaCa;(/) vanish. 

So for a given data point u, the cardinality of hber of tti over u , 7rj"^(u), measures 
the number of critical points. 

We want to discuss the locus of exceptional data points u at which the number 
of complex critical points is different from the ED degree. As mentioned in the 
introduction, there are three ways in which we can have different number of critical 
points then expected. The first one is because a critical point may wander off into 
SingA due to the closure appearing in the definition of Sx ■ This locus is called the 
ED data singular locus. 
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2.1. Data singular locus. We use the precise definition of the ED data singular 
locus from |3], that is 

7ri(£x n 7r^^(SingX)). 

We denote the ED data singular locus of an algebraic variety X by DS(X) (abbre¬ 
viating ’’data singular” locus) and we aim to describe the data singular locus of 
affine cones. Our main result in this section is the following theorem. 

Theorem 1. Let X C C" be an irreducible affine cone that is not a linear space. 
Then the following two inclusions hold 

X* DS(X) C( 2 ) X* + SingX, 

where X* denotes the dual variety to X. 

We view X* as subset of C” via the standard symmetric bilinear form (-j-) on 
C”. 


Proof. Eirst we prove inclusion (1) for a dense subset of X*. Eor this take u G X*, 
such that there exists a regular point Xr G X’'®®, such that u _L T^^X, that is all 

u 

J 3’Cxj. (d) 

X and JaCa;^(/) is the Jacobian of the (radical) ideal J of at the point Xr. We 

u 


the (c-l- 1) X {c+ 1) minors of 
X and JaCa;^(/) is the Jacobie 
denote an arbitrary (c J- 1) x (c -I- 1) minor of this matrix by 


vanish, where c is the codimension of 


J aCa;,. ( J) 


(c-ei) 


We claim that {u + Xxr, Xxr) G Sx for all real A > 0. We have that if / S J, 
homogeneous of degree d, then V/(Ax) = f{x). So if x^ is a regular point then 

\xr is also regular, for any A > 0. Moreover we get that for any (c -I- 1) x (c J- 1) 
minor 


f (w -I- Xxr) — XXr 
y JaC.\a,,. (/) 



JaCAx,(J) /(c+1) 


u 

J aCa;,, (d) 



= 0 , 


where N is the sum of degrees of the defining polynomials of I. 

So {u + Xxr, Xxr) G £x for all real A > 0. But then taking the limit when A goes 
to zero, we get that 

(u,0) G £xT TT^^SingX), 

since £x H 7r^^(SingX) is Zariski closed (hence closed wrt. Euclidean topology as 
well) and since 0 G SingX. Indeed, for every x G X the line {A • cc} is in the 
tangent space to 0, so TqX is equal to the the linear span of Jf, which has a grater 
dimension that X if and only if X is not a linear space, hence 0 G SingX. So then 
u = 7ri((u,0)) G DS(X). 

Eor the proof of (2) take an element (u, xq) G £x IX 7i^^(SingX), then this point 
can be approximated by a sequence in the part of £x over X’’®®. That is there 
exists a sequence —>■ 0 in C” and Xi -G xq with all the Xi G X"'®s, such that 


(u + Si,Xi) G £x- 

By the ED Duality Theorem for affine cones (see. |31 Theorem 5.2]) we get that 
(m -I- 6i) — Xi G X*, for all i. Now taking the limit, when i goes to infinity, we get 
that u — XqG X*, since X* is closed (hence closed wrt. Euclidean topology as well). 
Einally this means that u G Xq + X* C SingX + X*. □ 
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Note that the condition in the theorem that X is not a linear space is necessary 
in order to prove a similar statement. Otherwise if X is a linear subspace of C", 
then it has a non-empty dual (its orthogonal complement with respect to the inner 
product), but its singular locus is empty, hence its data singular locus is empty as 
well. 

2.2. Data isotropic locus. A second possibility for a data point u to have smaller 
number of critical points than expected is by letting one of the critical points to 
become isotropic. The authors of |3] define the ED degree of a projective variety 
in to be the ED degree of the corresponding affine cone in C”, moreover 

given a data point u the critical points to these two objects are in a one-to-one 
correspondence, given that non of the critical points lies in the isotropic quadric 
(see[31 Lemma 2.8]). In particular, the role of Q exhibits that the computation 
of ED degree is a metric problem. This is the reason that even though in the 
definition of the affine Ex we keep the isotropic critical points, but when we pass 
to projective varieties we will exclude the isotropic points. This way the data 
isotropic locus represents the locus of data points which have different number of 
critical points if X is considered as an affine cone or if is considered as a projective 
variety. 

More precisely we define the ED data isotropic locus to be 

tti{Ex niT^^iQnx), 

where Q = (X]r=i denotes the isotropic quadric with respect to the standard 
symmetric bilinear form. 

We denote the ED data isotropic locus of an algebraic variety X by DI(X) 
(abbreviating ’’data isotropic” locus). We have the following theorem for the ED 
data isotropic locus of affine cones. 

Theorem 2. Let X C C" be an irreducible affine cone. Then the following two 
inclusions hold 

X*C( 1 ) di(x)C( 2 ) x* + iQnx), 

where X* denotes the dual variety to X. 

Again we view X* as subset of C" via the standard symmetric bilinear form (-j-) 
on C”. 

Proof. The proof follows the lines of the proof of[^ keeping in mind that 0 S X is 
always an isotropic point. □ 

In the following two sections we will give examples to show that both inclusions 
appearing in Theoremand Theorem]^ can be strict and/or equalities. 

3. Examples of the ED data singular locus 

In this section we present several useful examples concerning the ED data singu¬ 
lar locus of an affine cone. Before we get to the examples we present how can one 
computationally determine the objects we are working with. We illustrate the main 
algorithms with code in Macaulay2 [5]. For an affine cone X C C", of codimension 
c with defining radical ideal I, one can determine its dual X* using the following 
code by nn Algorithm 5.1]. 
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Example 3 (Computing the dual variety). We present the algorithm for the 
real affine cone X C defined by the homogeneous equation / = ccf + x^x^. 

n=3; 

kk=QQ[x_l..x_n,y_l..y_n] ; 
f=x_l"3+x_2~2*x_3; 

I=ideal(f); 
c=codim I; 

Y=matrix{{y_l..y_n}}; 

Jac= Jacobian gens I; 

S=submatrix(Jac,{0..n-l},{0..numgens(I)-l}); 

Jbar=SI transpose(Y); 

EX = I + minors(c+1,Jbar); 

SingX=I+minors(c,Jac); 

EXreg=saturate(EX,SingX); 

IDual=eliminate(toList(x_l..x_n),EXreg) 

Which gives at the end that X* is the zero locus of the polynomial f* = Axf— 27 X 2 X 3 

Following the dehnition of the data singular locus, the next example contains an 
algorithm for calculating the ideal of it. 

Example 4 (Computing the data singular locus). We present the algorithm 
for the real affine cone X C defined by the homogeneous equation / = xf + X 2 X 3 . 

n=3; 

kk=QQ[x_l..x_n,y_l..y_n] ; 
f=x_l"3+x_2~2*x_3; 

I=ideal(f); 
c=codim I; 

Y=matrix{{x_l..x_n}}-matrix{{y_l..y_n}}; 

Jac= Jacobian gens I; 

S=submatrix(Jac,{0..n-l},{0..numgens(I)-l}); 

Jbar=SI transpose(Y); 

EX = I + minors(c+1,Jbar); 

SingX=I+minors(c,Jac); 

EXreg=saturate(EX,SingX); 

DSX=radical eliminate(toList(x_l..x_n),EXreg+SingX) 

Which gives as output that DS(X) is the zero locus of the polynomial xi{'ix\ — 
27x\x3). 

Now we arrived at the point to present a sequence of interesting varieties and 
the corresponding duals and data singular loci. The first example is the one we 
used for presenting the algorithms previously. In this example both inclusion (1) 
and inclusion (2) are strict, as it will be seen. 

Example 5 (Cuspidal Cubic Cone). Let X C be the real variety defined by 
the homogeneous equation / = xf + X 2 X 3 . Since it is an affine cone it has a dual 
X*, which is dehned by the dual equation f* = Axf — 27 x 2 X 3 . For the data singular 
locus we get that DS(X) is the zero locus of the polynomial a;i(4a;f — 27 X 2 X 3 ). So we 
can see that X* is even a component of DS(X). Moreover X* +SingX is something 
much larger and not equal to DS(X). For example the point (3, 2,1) + (0, 0,1) € 
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X* + SingX, but is not on DS(X). Figure 1 shows X in blue and X* in green and 
DS(X) is the union of the green colored X* and the additional surface in red. 



Figure 1. V{x\ + X^x^) together with its dual and its data sin¬ 
gular locus 


The next example shows that both inclusions (1) and (2) can be in fact equalities. 
More generally we have the following corollary to Theorem 

Corollary 6 . Let X C C” he an ajfine cone, with SingX = {0}, then DS(X) = X*. 
Moreover if X is a general hypersurface of degree d, then 

deg(DS(X)) =d((i-l)”-i. 

Proof. The first part follows directly from the claim of Theorem The moreover 
part is classical and we refer to m Exercise 5.14]. □ 

Example 7 (Cone over ellipse). Let X C the cone over an ellipse, defined 
by the homogeneous equation f = xf + dxl — 9x§. The singular locus SingX only 
contains 0, so as a consequence of Theoremwe have that DS(X) equals the dual 
variety X*, defined by the dual equation f*=xl + a;|/4 — ai§/9. Figure 2 shows 
X in blue and X* in green. 

The next example concernes the well known and much used determinantal vari¬ 
eties. We will see that for this variety inclusion (1) is strict and inclusion (2) is an 
equality. 
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Figure 2. V { x 1 + 4x1 — 9x|) together with its dual 


Example 8 (Determinantal varieties). Denote by the variety of n x m 

matrices (suppose n < m) of rank at most r. It is classical that the singular locus 
is the variety . By |3 Chapter 1, Prop. 4.11] we have that the dual variety 

is exactly So applying Theorem]^ we get that 




) C 


So for rank-one matrices (r = 1) we get that DS(M^^^) = which is not 

a surprise based on Corollary since is smooth, except 0. But something 

more is true for general r. We claim that the upper bound for the inclusions is 
always attained. For this we have the following proposition. 


Proposition 9. The ED data singular locus of the determinantal variety is 

equal to , for alll <r <n — \ . 


Proof. A n X m matrix U lies in DS(M,y(^^) if and only if it has a singular critical 
point. By [31 Example 2.3] all the critical points of U look like 

Ti •Diag(0,0, ...,CTii,0, ...,0,(Ti,.,0, ...,0) ■T 2 , 

where the singular value decomposition of U is equal toU = Ti •Diag((Ti,..., <Jn)-T 2 , 
with a I > ... > CTn singular values and Ti,T 2 orthogonal matrices of size nxn and 



8 


EMIL HOROBET 


m X m. Such a critical point is singular if and only if it has rank at most r — 1, 
which can only happen if one of the singular values is zero. So there 

exists a singular critical point to U if and only if there is a zero singular value of 
U, which can only happen if U has a rank defect, hence all the (n — 1) x (n — 1) 
minors are zero, that is 17 G □ 

The next example shows that X* is a subvariety of DS(7f) but not necessarily a 
component of it. 

Example 10 (Hurwitz determinant). In control theory, to check whether a 
given polynomial is stable one builds up the so called Hurwitz matrix and 
checks if every leading principal minor of is positive. Take n = 4, then the 4-th 
Hurwitz matrix looks like 


/ 

2:2 

2:4 

0 

0 

\ 


2:1 

2:3 

2:5 

0 



0 

2:2 

2:4 

0 


V 

0 

2:1 

2:3 

2:5 

/ 


The ratio r 4 = det(i 74 )/a ;5 is a homogeneous polynomial and it is called the Hurwitz 
determinant for n = 4 by [31 Example 3.5]. 

Let 77 C C® be the affine cone defined by r 4 . Then its dual variety has one 
irreducible component given by 

X* = V^—xzXi -\- X 2 Xz, —xl -I- X 1 X 5 , —X 2 X 3 -I- a;ia;4). 

While its data singular locus DS(X) has two irreducible components and it is defined 
by 

V{{xix'^ + X2X^Xi -\- x'Ix 5){X2X3 — X1X2X4 — 2XiX2X^ — X3X4 + 2X2X4X5 + 0:23^42:5)). 

It is clear that X* is not component of DS(77). Moreover DS(X) is not equal to 
X* + SingX, since Sing77 = V{x 2 ,X 4 ) and the point 

(2,1,1,0,1) = (1,1,0,0,0) + (l,0,l,0,l) 

lies on X* + Sing77 but it is not on DS(77). 

We have thus seen examples of varieties with: both inclusions in Theorem[2being 
strict, both inclusions in Theorembeing equalities and the second inclusion being 
an equality, while the first one is strict. It is natural to ask if there are examples 
where the first inclusion is an equality, while the second one is strict. The author 
could not find such an example, so the following question arises. 

Problem 11. Find an affine cone X, such that X* = DS C 77*-|-Sing(77) or prove 
that there is no such X. 

4. Examples of the ED data isotropic locus 

In this section we present several application oriented examples concerning the 
ED data isotropic locus of an affine cone. We begin with presenting how can one 
computationally determine the data isotropic locus of a variety. 

Example 12 (Computing the data isotropic locus). We present the algorithm 
for the real affine cone 77 C C® defined by the homogeneous equation / = XiXq — 
X2X3 -b 2 : 32 : 4 , representing the Grassmanian of planes in 4-space. 




THE DATA SINGULAR AND THE DATA ISOTROPIC LOCI FOR AFFINE CONES 


9 


n=6; 

kk=QQ[x_l..x_n,y_l..y_n]; 
f=x_l*x_6-x_2*x_5+x_3*x_4; 

I=ideal(f); 
c=codim I; 

Y=matrix{{x_l..x_n}}-matrix{{y_l..y_n}}; 

Jac= Jacobian gens I; 

S=submatrix(Jac,{0..n-l},{0..numgens(I)-l}); 

Jbar=SI transpose(Y); 

EX = I + minors(c+1,Jbar); 

SingX=I+minors(c,Jac); 

q=sum for i from 1 to n list x_i"2; 

Q=ideal(q); 

EXreg=saturate(EX,SingX); 

DIX=radical eliminate(toList(x_l..x_n),EXreg+Q) 

Which gives that DI(X) is the zero locus of the polynomial XiXq — X 2 X 5 + X 3 X 4 , so 
we get that the data isotropic locus is equal to the dual variety which equals the 
variety. 

The next example shows that the data isotropic locus can be equal to the dual 
and strictly contained in X* + {X D Q). 

Example 13 (Cayley-Menger variety). Let X denote the variety in with 
parametric representation 



among 3 interacting agents with coordinates zi, Z 2 and z^ on the line K. The prime 
ideal of X is given by the determinant of the Cayley-Menger matrix 



So X is defined by the irreducible polynomial 

f = x\ — 2 XiX2 x\ — 2 X 1 X 3 — 2 X 2 X 3 + X3. 


After running the computations one can see that the data isotropic locus equals 
the dual variety, which is defined by /* = X 1 X 2 + X 1 X 3 + X 2 X 3 , see Figure 3. And 
it does not equal X* {Q Cl X) for example because the point (1, 0, 0) + (0,1, i) G 
X* + (Q n X), but it does not lie on DI(X). 


The next example shows that both inclusions from Theoremcan be strict. 


Example 14 (Cayley’s Cubic). Let X be defined hy f = x\ — xix\ — X 1 X 3 + 
2 x 2 X 3 Xi — the 3x3 symmetric determinant in C^. This hypersurface is 

sometimes called the Cayley’s cubic surface and receives much attention in the 
study of elliptopes and exponential varieties in algebraic statistics, see for instance 
[TTl Example 5.44], (TUI Example 1.1] and [5]. Its dual variety is the quartic Steiner 
surface defined by f* = cclxg — 2a:ia;2a;3X4 + x^x^ + XgX^. After running the 
computations one finds that the data isotropic locus is the union 


DI(Ar) = E(x]® + 4x]®X2 + 6 x 



... + 12 ^x\x\^)\jX*. 
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Figure 3. Cayley-Menger variety (in blue) together with its dual 
(in green). 


So it is clearly not equal to the dual variety. And it is not equal to the X* + {QHX) 
either, because for example the point 

(i,i,0,0) + (0,0,i,i) e A* + (QnA), 

but it is not in DI(A). 

Our next example shows that the second inclusion in Theorem can be equality 
and moreover it can give the whole space. 


Example 15 (Special essential variety). Essential matrices play an important 
role in multiview geometry, see for instance [7]. The connections between the ED 
degree theory and multiview geometry were investigated in [3l Example 3.3]. The 
set of essential matrices is called the essential variety and it is defined as follows 

£ = {X £ Maxal det A = 0, 2XX^X - trace(AA^)A = 0}. 


It is a codimension 3 variety of degree 10. The ED degree of £ is 6, as was proved 
in [H Example 5.8]. We are interested in the data isotropic locus of this variety, 
but because of computational reasons we will take a linear section of it and we will 
only consider the symmetric, constant diagonal essential matrices , which we will 
call the special essential variety and will denote by S£. More precisely we define 
S£ to be 


A = 


Xi 

X2 

X 3 

X2 

Xi 

Xi 

X3 

X 4 

Xi 


det A = 0, 2AA^A - trace(AA^)A = 0 
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Since this variety will not be irreducible we will carry out our computations com¬ 
ponentwise. When running the computations one will find that the data isotropic 
locus is the hole space. Indeed one can observe that SS is inside the isotropic 
quadric Q, so every critical point is isotropic. We have that 

DI(X) =X* + {Xf^Q)=X* +X = <[:'^. 

Moreover DI(X) is not equal to the dual variety, since X* is a proper variety dehned 
by f* = {x 2 + x‘l){x 2 + x‘^){xl + x^. Moreover it is clear that the dual is not a 
component of DI(X). 

In the last example the reader can see that both inclusions from Theorem can 
be equalities. 

Example 16 (Line through the origin). In what follows let X be a line through 
the origin in C^. So we have X = V{xi + 2 x 2 + 3a;3,4a:i -I- 5x2 + 6x3). Then we 
get that X intersects the quadric Q only in the point 0, so by Theorem we 
get immediately, that X* = DI(X) = X* + {0}, and the dual is the orthogonal 
complement of X, so it is defined by xi — 2x2 + X 3 . 
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